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Given a set of independent Poisson random vari- 
ableswith common mean, we study the distribution of 
their maximum and obtain an accurate asymptotic for- 
mula to locate the most probable value of the maximum. 
We verify our analytic results with very precise numer- 
ical computations. 



We deal with a set of independent Poisson random variables {Xi, X2, . . . , Xn} 
with common mean A, so that Pr[Xj — k] — e'^X^ /k\. We let M„ = max(Xi) 
and wish to describe the distribution of M^. Our motivation is a problem 
in random graph theory, where we were interested in the distribution of 
maximum degree in graphs with Poisson degree distribution. 
We have 



where Q and r(-, •) are incomplete Gamma functions; that is. 



Pr[M„ ^k]^ Pr[Xi ^ A;]" - Q(A;+1, A)" = r(fc+l, A)"/r(A;+l)". 



Our aim is to approximate the distribution of Mn- We have 



Pr[Xi < A;] = Q{k,\) = r{k,X)/r{k) 




Prom the independence of the Poisson variables. 



Pr[M„ = k] 



Pr[M„ ^ A;]-Pr[M„ ^ 
Q(A;+l,Ar-g(A;,Ar 
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Figure 1: From left to right: the distribution of the maximum of Poisson 
variables for A = 1/2,1,2,5 (left to right) and n = 10°, 10^, 10^ . . . , lO^^. 
Note that there is an error in Fig 1 in Anderson et al. (1997), where the 
curves labelled k = 6 and k = 8 are incorrect. 
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Figure 2: The maximal probability (with respect to /.„) that M„ G {/„, 

for A = 1/2, 1,2,5 (left to right) and 10° ^ "n. ^ 10^°. The curves show the 

probability that M„ takes either of its two most frequently occurring values. 

Examples of these distributions are shown in Figure 1. These numerical 
results demonstrate the so-called focussing effect; the maxima M„ are con- 
centrated on at most two adjacent integers for large n; we call them modal 
values. It is this focussing that allows us to characterize the distributions 
very precisely by a single asymptotic estimate. 

In previous work on this problem, Anderson (1970) proved the existence of 
integers /„ such that Pr[M„ G (/„, /n+1)] — > 1 as n — > oo for fixed A > 0; and 
that In ~ /3n, where is defined as the unique solution of Q{Pn, A) = l/n. 

Following this work, Kimber (1983) computed an asymptotic result; he 
showed /„ ~ logrz/loglogn and P„ ~ {k/InY^^'' with_B„ dense in [—1/2, 1/2], 
and concluded that to the first order, the rate of growth of is inde- 
pendent of the Poisson parameter A. He concluded that P„, defined as 
Pn = Pr[M„ G (/n, /n + 1)], osciUatcs and the oscillation persists for arbitrar- 
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ily large n. We illustrate in Figure 2 exactly how this probability oscillates 
as n ^ oo. Our numerical experiments show that log n/ log log n estimates 
In very poorly. We aim to improve this asymptotic formula. 

Our method is a refinement of that of Kimber; that is, we consider a 
continuous distribution g which interpolates the Poisson maximum distribu- 
tion, and we solve g{x) = 1/n. Consider g\{x) = 1 — r(x + l, A)/r(x + l) for 
fixed A G M"*", which is a strictly decreasing function on (0, oo). If e = 1/n 
is a small positive real, then g\{x) has a unique root x(e) which increases as 
e ^ 0"*". We will develop an asymptotic expansion (as e — 0) of this root 
x(e). 

We have 

°° A* 

gx{x) = exp(-A) A^ 



T{x + t + l] 



and we will work with 
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\og{gx{x)) = -a;log(a;) + (l + logA)a;--log(a;) 

A first approximation to the solution of log{gx{x)) = — logra large and neg- 
ative is given by keeping only the dominant first two terms in Equation (1): 

M„ ~ Xq 



f logn \ 
\^exp(l)A ) 



w 



where W{-) is the principal branch of Lambert's W function (Corless et al., 
1996). That this is already quite accurate can be seen form the dark blue 
curves in Figure 3. However, we would like to do better; ideally the error 
should be less than unity so that the mode of the distribution is correctly 
identified. A refinement xi may be generated by making a single Newton 
correction step; that is, xi = xo — {h{xo) + \ogn)/h'{xo), where h is some 
approximation to \og{gx)- For example, by keeping all terms in log{gx{x)) 
and \og{gx{x)y which do not vanish as n — * oo, we obtain 

logA-A-log(27r)/2-31og(xo)/2 

Mn^ Xi= XoH . r . . 

log(a;o)-logA 

This appears to have error less than unity for all values of n and A considered 
in Figure 3, and so is probably sufficiently precise for all practical purposes. If 
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Figure 3: Exact values and asymptotics of /„ for A = 1/2,1,2,5 (left to 
right) and n = 10°, . . . , 10^°. The staircase red line (almost hidden by the 
cyan line) represents the exact mode /„; the other hues represent asymptotic 
approximations: green for the result of Kimber (1983) (which is independent 
of A), dark blue and cyan for our new results xq and xi respectively. The 
cyan curve always sits between the steps of /„, meaning that xi has error 
less than unity. 

further accuracy is needed, it may be obtained by additional Newton steps. 
In any case, both xq and xi are considerably more precise than Kimber's 
approximation. 
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